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Abstract 

A measurement of the cross-section for pp — t- Z — t- e'''e~ is presented using data at 
y/s = 7TeV corresponding to an integrated luminosity of 0.94 fb~^. The process is 
measured within the kinematic acceptance pT > 20 GeV/c and 2 < ry < 4.5 for the 
daughter electrons and dielectron invariant mass in the range 60-120 GeV/c^. The 
cross-section is determined to be 

cj(pp ^ Z ^ e+e") = 76.0 ±0.8 ±2.0 ±2.6 pb 

where the first uncertainty is statistical, the second is systematic and the third is 
the uncertainty in the luminosity. The measurement is performed as a function of 
Z rapidity and as a function of an angular variable which is closely related to the 
Z transverse momentum. The results are compared with previous LHCb measure- 
ments and with theoretical predictions from QCD. 
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1 Introduction 



The measurement of vector boson production permits a number of tests of electroweak 
physics and of quantum chromodynamics (QCD) to be performed. In particular, the 
angular acceptance of LHCb, roughly 2 < rj < 5 in the case of the main tracking sys- 
tem where rj denotes pseudorapidity, complements that of the general purpose detectors 
ATLAS and CMS. LHCb measurements provide sensitivity to the proton structure func- 
tions at very low Bjorken x values where the parton distribution functions (PDFs) are 
not particularly well constrained by previous data from HERA (see for example Ref. [1]). 

The most straightforward decay modes in which the and Z bosons can be studied 
using the LHCb data are the muonic channels, Z — )■ and W"*" — )• \L^y^- Measure- 

ments of Z — )■ (J-^M-^ and of Z — ^ t+t^ using the LHCb data at -y/s = 7TeV have already 
been presented [2 |3]. To complement these studies, the electron channels Z — e"'"e~ 
and W"*" — )■ e+Ve, which offer statistically independent samples with different sources of 
systematic uncertainties, are examined. 

The main difficulty with electroiij^ reconstruction in LHCb is the energy measurement. 
A significant amount of material is traversed by the electrons before they reach the mo- 
mentum analysing magnet, and their measured momenta are therefore liable to be reduced 
by bremsstrahlung. For low energy electrons, the bremsstrahlung photons can frequently 
be identified in the electromagnetic calorimeter and their energies added to the measured 
momentum of the electron. However, in the case of and Z decays, the electrons are of 
high momentum and transverse momentum (pt), so that the bremsstrahlung photons of- 
ten overlap with the electrons. The LHCb calorimeters were designed so as to optimise the 
the measurement of photons and 7t°s from heavy fiavour decays, whose transverse energy 
[Et) values are generally well below 10 GeV. As a consequence, individual calorimeter 
cells saturate at around 10 GeV, so it is not possible to substitute the calorimeter en- 
ergy for the momentum measured using the spectrometer. We therefore have a situation 
in which the electron directions are well determined, but their energies are underestimated 
by a variable amount, typically around 25%. Nevertheless, the available information can 
be used to study certain interesting variables. 

In this paper, we present a measurement of the cross-section for pp — j- Z — )■ e"'"e~ 
using the data recorded by LHCb in 2011 at y/s = 7TeV. Throughout this paper we 
use Z — )• e^e~ to refer to the process Z/7* — i- e+e^ where either a virtual photon or a 
Z boson is produced and decays to e+e~. For consistency, the measurement is presented 
in the same kinematic region as the recent measurement of Z — )■ (J-'''IJ.~ using the 2010 
LHCb data at y/s = 7TeV [2j: 2 < < 4.5 and px > 20 GeV/c for the leptons and 
60 < M < 120 GeV/c^ for the dileptons where M is the invariant mass. Since the rapidity 
of the Z boson can be determined to a precision of ~0.05, the rapidity distribution will 
be presented. However, the of the Z boson is poorly determined and its distribution 
will not be discussed. A similar problem was encountered by the DO collaboration 

"'^The term "electron" is used generically to refer to either e"*" or e~. 
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who employed a new variable proposed in Ref. depending only on track angles 



where M and refer to the lepton pair, A?7 and A0 are the differences in pseudorapidity 
and azimuthal angles respectively between the leptons, and the acoplanarity angle is 
0acop = VT — |A0|. The pt of the Z boson is correlated with 0*, and the resolution on 0* 
is excellent, with a precision better than 0.001. The measurement of 0* presented here 
therefore largely accesses the same physics as a measurement of the Z distribution. The 
measurement of the distribution of Z rapidity (denoted yz) is expected to show sensitivity 
to the choice of PDFs, while 0* is likely to be more sensitive to higher order effects in the 
QCD modelling. 

After a brief description of the detector. Sect. |3] describes the event selection, and 
Sect. |4] outlines the determination of the cross-section. The results are given in Sect. [5] 
followed by a short summary. 



The LHCb detector [6] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < 77 < 5, designed primarily for the study of particles containing b or c quarks. 
The detector includes a high precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4Tm, and three stations 
of silicon-strip detectors and straw drift tubes placed downstream. The combined track- 
ing system has a momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% 
at lOOGeV/c for hadrons and muons, and an impact parameter resolution of 20 |am for 
tracks with high transverse momentum. Charged hadrons are identified using two ring- 
imaging Cherenkov detectors. Photon, electron and hadron candidates are identified by 
a calorimeter system consisting of scintillating-pad (SPD) and preshower (PRS) detec- 
tors, an electromagnetic calorimeter (ECAL) and a hadronic calorimeter (HCAL). The 
acceptance of the calorimeter system is roughly 1.8 < 77 < 4.3. Muons are identified by a 
system composed of alternating layers of iron and multiwire proportional chambers. 

The trigger [7j consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage which applies full event reconstruction. A 
significant improvement to the trigger was implemented during August 2011 which affected 
the trigger efficiency for Z — )■ e"'"e~. The data samples before and after this change are 
treated separately and will be referred to as data sample I and data sample II. These 
correspond to integrated luminosities of 581 ± 20 pb^^ and 364 ± 13pb~^ respectively, 
yielding a total of 945 ± 33pb~^. 
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2 LHCb detector 
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3 Event selection 



The Z — )■ e"'"e~ sample is initially selected by single-electron triggers, which require elec- 
trons to have an Et above a given threshold between 10 and 15 GeV depending on the 
data-taking period and specific trigger. The Z — )■ e"'"e~ selection starts from a sample 
of e"'"e~ candidates with high invariant mass, which is refined by requiring the following 
selection criteria: 

• At least one of the candidate electrons must be selected by a high-i^x electron 
trigger. 

• The electrons are both required to have pt > 20 GeV/ c and pseudorapidity in the 
range 2.0 < rj < 4.5. The invariant mass of the e"'"e~ pair should be greater than 
40 GeV/c2. 

• Requirements on calorimeter information are imposed to provide particle identifica- 
tion (PID) of electrons. The particle must satisfy Eecal/pc > 0.1, where p is the 
particle momentum, with bremsstrahlung correction if available, and ^'ecal is the 
ECAL energy associated with the particle. The particle is required to lie within the 
HCAL acceptance and to satisfy Eucal/pc < 0.05, where -Ehcal is the HCAL en- 
ergy associated with the particle. The energy in the preshower detector associated 
with the particle is required to satisfy -EpRs > 50 MeV. These requirements im- 
pose an electromagnetic shower profile, while being loose enough to maintain a high 
electron efficiency despite the effects of calorimeter saturation and bremsstrahlung. 

• If more than one Z — )■ e"^e~ candidate satisfies the above requirements in an event, 
just one candidate is used, chosen at random. This only affects around 0.5% of 
cases, and in all instances the multiple candidates share one daughter. 

A sample of same-sign e^e^ combinations, subject to the same selection criteria, is used 
to provide a data-based estimate of background. The main background is expected to 
arise from hadrons that shower early in the ECAL and consequently fake the signature of 
an electron. These will contribute approximately equally to same-sign and opposite-sign 
pairs. The contribution from semileptonic heavy flavour decays should be similar to the 
small level (~ 0.2%) estimated for the Z — )■ channel [2]; in any case, subtracting 

the same- sign contribution should account for most of this effect. 

Simulated event samples of Z — )■ e"'"e~ with M(e"'"e~) > 40 GeV/ are also used to 
assess some efficiencies as discussed below. Simulated samples of Z — )■ t+t~ and of tt are 
used to assess possible background contributions. For the simulation, pp collisions are 
generated using Pythia 6.4 [8j with a specific LHCb configuration [9J and the CTEQ6L1 
PDF set [in] • The interaction of the generated particles with the detector and its response 
are implemented using the Geant4 toolkit [H] as described in Ref. [T2]. Simulated 
samples based on different versions of GEANT and of the detector model are employed, 
which allows the reliability of the simulation to be assessed. The simulated events are then 
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Figure 1: Invariant mass distribution of Z — )■ e"*"e~ candidates. The data are shown as 
points with error bars, the background obtained from same-sign data is shown in red 
(dark shading), to which the expectation from signal simulation is added in yellow (light 
shading). The Z — )• e"'"e~ simulated distribution has been normalised to the (background- 
subtracted) data. 



reconstructed in the same way as the data, including simulation of the relevant trigger 
conditions. 

The invariant mass distribution of the selected candidates is shown in Fig. [1} The 
distribution falls off abruptly above the Z mass and is spread to lower masses by 
bremsstrahlung. Good agreement in shape is observed between data and the simulation 
sample used in the data correction; this will be further discussed below. The background 
estimated from same-sign events amounts to 4.5% of the total number of e"'"e~ candi- 
dates. The backgrounds from t+t" and tt events are estimated to be around 0.1% and 
are neglected. 



4 Cross-section determination 

In a given bin of Z rapidity or 0*, the cross-section is calculated using 

I V + -^ N(e+e-) - N(e±e±) „ „ 

(t(pp ^ Z ^ e+e ) = ^ '- ^ ■ fpsR ■ fMZ , (2) 

CGEC ■ ^trig ■ Ctrack ' ^kin ' ^PID ' J ''-Clt 

where iV(e"'"e~) is the number of Z candidates selected in data, A^(e^e^) is the background 
estimated from the number of same-sign candidates and J £dt is the integrated luminosity. 
The cross-section cr(pp — )■ Z — )■ e"'"e~) denotes the product of the inclusive production 
cross-section for the Z or 7* and the branching ratio to e"'"e~ . The meaning and estimation 
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Table 1: Quantities entering into the cross-section determination, averaged over the range 
of Z rapidity used. 





Data sample I Data sample II 




581 ±20 364 ±13 


cgec 


0.947 ±0.004 


Ctrig 


0.715 ±0.021 0.899 ±0.003 


Ctrack 


0.913 ±0.015 


Ckin 


0.500 ±0.007 


CPID 


0.844 ±0.011 


/fsr 


1.049 ±0.005 


/mz 


0.967 ±0.001 



of the other factors are described below. The values obtained for each, averaged over the 
acceptance, are summarised in Table [T] 

The luminosity is determined as described in Ref. [13] and has an uncertainty of 3.5%. 
The factor /fsr accounts for the effects of final-state electromagnetic radiation, correcting 
the measurement to the Born level. As in the Z ^ [i^\Ji~ analysis |2j it is determined 
using Photos [H] interfaced to Pythia [S], with Horace [T3] used as a cross-check. An 
overall systematic uncertainty of 0.5% is assigned to this correction [I6]. The factor /mz 
corrects for e"'"e~ events outside the mass range 60 < M(e^e~) < 120 GeV which pass the 
event selection, and is estimated from simulation by examining the true mass for selected 
events. 

The probability for a Z — )• e"'"e~ event to satisfy the trigger and selection requirements 
is given by the product of the efficiency factors, e, as described below. 

• Global event cuts (GEC) are applied in the trigger in order to prevent very large 
events from dominating the processing time. Their efficiency for selecting signal 
events is given by ecEC- In the Z — )■ e"'"e~ case, the most important requirement is 
on the multiplicity of SPD hits, A^spd < 600. This is strongly correlated with the 
number of primary vertices reconstructed in the event. The inefficiency is assessed 
by comparing with Z — )■ fJ-'^M-" candidates recorded in the same running period using 
a dimuon trigger for which a less stringent requirement of 900 hits is imposed. A 
correction is made for the small difference in the numbers of SPD hits associated with 
the electrons and muons themselves. This procedure is adopted for each number of 
reconstructed primary vertices and the results are combined to obtain the overall 
efficiency. 

• The trigger efficiency for events passing the final selection, etrig, is determined from 
data. A sample of events triggered independently of the e"*" is identified and used 
to determine the efficiency for triggering the e"*", and likewise for the e~. Using 
the total numbers of candidates for which the single electron trigger is satisfied at 
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each stage by the e"*" (A^"^), by the e~ (A^~) and by both {N^ ), the efficiency for 
triggering the e"*" is given by e"*" = /N~. The overall efficiency is then taken 
to be + — e^e^ assuming that the e+ and e~ are triggered independently. 
The procedure is validated on simulated events. The determination is performed 
separately in each bin of Z rapidity and 0*. In all cases, the statistical uncertainty 
on the efficiency is taken as a contribution to the systematic uncertainty on the 
measurement. 

• The track-finding efficiency, etrack, represents the probability that both of the elec- 
trons are successfully reconstructed. The simulation is used to determine the track- 
finding efficiency, in bins of Z rapidity and 0*, by calculating the probability that, 
in a Z — )■ e"'"e~ event whose generated electrons lie within the kinematic acceptance, 
both of the electrons are associated with reconstructed tracks that satisfy the track 
quality requirements, but not necessarily the kinematic requirements. Its statistical 
precision is propagated as a contribution to the systematic uncertainty. 

This efficiency is checked in data using a tag-and-probe approach. One electron is 
tagged using the standard requirements, and a search is made for an accompanying 
cluster of electromagnetic energy having a high and forming a high invariant 
mass with the tag electron. If such a cluster has no associated track it provides 
evidence of a failure to reconstruct the other electron. This sample contains signif- 
icant background, which can be discriminated by examining the distribution of 
the tag electron for cases where the photon candidate is and is not isolated. The 
Pt distribution of the electrons in signal events in data displays a clear shoulder 
extending to ~ 45 GeV/ c while that for background falls monotonically, as shown in 
Fig. [2j The number of signal-like events in which a cluster is not associated with 
a track can be used to estimate a tracking efficiency, and the ratio of efficiencies 
between data and simulation is applied as a correction to the tracking efficiency. 
The precision of the test is taken to define a systematic uncertainty, assumed to be 
fully correlated between bins of rapidity and 0*. 

• The kinematic efficiency, ekin, represents the probability that, in a Z — )■ e^e~ event 
whose generated electrons lie within the kinematic acceptance and are associated 
with reconstructed tracks, both tracks pass the kinematic selection requirements 
2 < f] < 4.5 and pt > 20GeV/c. The efficiency is estimated from simulation, with 
its statistical precision being treated as a contribution to the systematic uncertainty. 
This determination relies on a correct simulation, which can be tested using data. 
For example, underestimation of the amount of material in the simulation would 
cause a discrepancy between data and simulation in the pt distributions of the 
electrons or the reconstructed mass spectrum shown in Fig. [T] By comparing the 
shapes of the reconstructed mass spectrum and other kinematic distributions in 
data with different simulation samples, a systematic uncertainty on the momentum 
scale and hence on the kinematic efficiency is assigned. This is combined with the 
statistical uncertainty mentioned above, with the systematic contribution taken to 
be fully correlated between bins of rapidity and 0*. 
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20 30 40 50 60 



[GeV/c] 

Figure 2: Distribution of for the "tag" electron in cases where an isolated cluster 
of energy of high E^^ is found in the electromagnetic calorimeter. This is fitted with 
two components obtained from data, the Z — )■ e"'"e~ signal whose shape is taken from 
those candidates where the cluster is associated with an identified electron track, and 
background whose shape is obtained from candidates where the cluster is not isolated. 

• The PID efficiency, epiD, represents the probability that, in a Z — )■ e"''e~ event with 
reconstructed electron tracks satisfying the kinematic requirements, both tracks ful- 
fil the calorimeter energy requirements for identified electrons. This includes the 
probability that the tracks are within the calorimeter acceptance and have been 
successfully associated with calorimeter information. Because of the acceptance 
contribution, the efficiency has a strong dependence on the Z rapidity. This depen- 
dence is taken from simulation, while the overall normalisation of the PID efficiency 
is estimated directly from data, using a tag-and-probe method. 
Starting from a sample which requires just one high electron, events are selected 
by applying the usual criteria except that only one of the e"*" and e~ (the "tag") 
is required to pass the calorimeter-based electron identification requirements. The 
other track is used as a "probe" to test the PID efficiency. The requirement of only 
one identified electron admits a significant level of background, which is assessed 
similarly to the tracking efficiency by examining the distribution of the tag or 
alternatively the px of the probe electron or the invariant mass of the two particles. 
The size of the signal component can be used to define the number of Z events which 
fail the PID, and hence to determine the PID efficiency and its uncertainty. 

A systematic uncertainty is also assigned to the same-sign background subtraction. 
The assumption that same-sign e^e^ combinations model background in e^e~ events is 
tested by selecting events which satisfy all criteria except that one of the particles fails 
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the calorimeter energy requirements. This sample should be dominated by background, 
and shows an excess of ~8% of opposite-sign events over same-sign events. Accordingly 
a systematic uncertainty amounting to 8% of the number of same-sign events is assigned 
to the measurements. 



5 Results 

Using the efficiencies described above, the event yields detailed in Table |2] and Eq. ^ 
separate cross-section measurements for the two data-taking periods are obtained. Since 
these are in good agreement, the results are combined using a weighted average, and as- 
suming their uncertainties are fully correlated apart from the statistical contribution and 
the uncertainty in the trigger efficiency. Data sample II has a smaller integrated lumi- 
nosity but a higher and more precisely estimated trigger efficiency. The weighting of the 
two samples is chosen to minimise the total uncertainty on the cross-section integrated 
over Z rapidity. The values of the differential cross-sections obtained are given in Table |2] 
Correlation matrices may be found in the Appendix. The bin 4.25 < yz < 4.5 is empty in 
data, and is expected to have close to zero detection efficiency since the calorimeter ac- 
ceptance extends only slightly beyond 4.25. Hence no measurement is possible. However, 
the QCD calculations discussed below predict a cross-section below ~0.01pb in this bin, 
which is negligibly small, so comparisons with the Z— )■ (J-^M-^ results or with theoretical 
calculations in the range 2 < ?/z < 4.5 are still meaningful. 

The cross-section integrated over Z rapidity is obtained by summing the cross-sections 
of all bins of yz, taking the uncertainties associated with the GEC and the luminosity 
to be fully correlated between bins, along with parts of the tracking, kinematic and PID 
efficiencies, and treating the other contributions as uncorrelated. The cross-section is 
measured to be 

(t(pp Z ^ e+c") = 76.0 ± 0.8 (stat.) ± 2.0 (syst.) ± 2.6 (lumi.) ± 0.4 (FSR) pb, 

where the first uncertainty is statistical, the second is the experimental systematic un- 
certainty, the third is the luminosity uncertainty and the last represents the uncertainty 
in the FSR correction. Since the results have been corrected to the Born level using the 
factor /fsr, it is possible to compare this measurement with that found in the Z— )■ |J."^|J.~ 
analysis f2\ using 37pb~"^ of data, namely 76.7 ± 1.7 (stat.) ± 3.3 (syst.) ± 2.7 (lumi.) pb. 
Accounting for correlated uncertainties, the ratio of cross-sections is 

(t(pP Z — )■ e+e~) , , , , 

' = 0.990 ±0.024 (stat.) ±0.044 (syst.). 



cr(pp — )• Z — )■ |J,+ |J,~ 

This may be regarded as a cross-check of the analyses. Assuming lepton universality, 
the two cross-sections can be combined in a weighted average so as to minimise the total 
uncertainty, yielding 

(t(pP ^ Z ^ tr ) = 76.1 ± 0.7 (stat.) ± 1.8 (syst.) ± 2.7 (lumi.) ± 0.4 (FSR) pb. 
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Table 2: Event yields and measurements for the differential cross-section of 
pp — )■ Z — )■ e"'"e~ at a/s =7 TeV as a function of Z rapidity, |/z, and of 4>*. The first 
uncertainty is statistical, the second and third are the uncorrelated and correlated ex- 
perimental systematic uncertainties respectively, and the fourth is the uncertainty in the 
FSR correction. The common luminosity uncertainty of 3.5% is not explicitly included 
here. The results are given for the combined data sample. The right-hand column gives 
the values used for the FSR correction factor. 



yz 




e+e") 




±e±) 


da/dyz [pb] 


/fsr 


2 


00- 


-2 


25 




988 




40 




13.6 ±0.7 ±0.4 ±0.3 ±0.1 


1 


049 ±0 


004 


2 


25- 


-2 


50 




3064 




121 




39.4 ± 1.0 ±0.6 ±0.8 ±0.2 


1 


046 ±0 


002 


2 


50- 


-2 


75 




4582 




202 




56.7 ± 1.2 ±0.7 ± 1.3 ±0.3 


1 


050 ±0 


002 


2 


75- 


-3 


00 




5076 




214 




63. 2 ±1.3 ±0.8 ±1.5 ±0.3 


1 


049 ±0 


002 


3 


00- 


-3 


25 




4223 




181 




59.9 ±1.4 ±0.8 ±1.6 ±0.3 


1 


056 ±0 


002 


3 


25- 


-3 


50 




2429 




135 




43.8 ± 1.3 ±0.8 ± 1.1 ±0.2 


1 


054 ±0 


003 


3 


50- 


-3 


75 




906 




61 




20.5± 1.0±0.7±0.6±0.1 


1 


030 ±0 


006 


3 


75- 


-4 


00 




143 




18 




5.9 ±0.8 ±0.5 ±0.3 ±0.1 


1 


074 ±0 


029 


4 


00- 


-4 


25 




9 




2 


0.66 ± 0.44 ± 0.30 ± 0.04 ± 0.02 


1 


074 ±0 


029 


4 


25- 


-4 


50 
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* 






e+e") 


N{e 


±e±) 




da/d(f)* [pb] 




/fsr 







00- 


-0 


05 




9696 




363 




693 ±10±6±17±3 


1 


059 ±0 


001 





05- 


-0 


10 




4787 




219 




326±7±4±8±2 


1 


047 ±0 


002 





10- 


-0 


15 




2382 




115 




164 ± 5 ± 3 ± 4 ± 1 


1 


039 ±0 


002 





15- 


-0 


20 




1384 




80 




99.1 ±4.0 ±2.0 ±2.2 ±0.5 


1 


043 ±0 


003 





20- 


-0 


30 




1434 




82 




49.6 ±2.0 ±1.1 ±1.0 ±0.3 


1 


042 ±0 


003 





30- 


-0 


40 




707 




39 




25.5 ± 1.4 ±0.8 ±0.6 ±0.1 


1 


049 ±0 


004 





40- 


-0 


60 




583 




41 




10.8 ±0.7 ±0.4 ±0.3 ±0.1 


1 


052 ±0 


005 





60- 


-0 


80 




217 




13 


4 


05 ±0.38 ±0.20 ±0.09 ±0.03 


1 


054 ±0 


005 





80- 


-1 


00 




91 




9 


1 


41 ±0.23 ±0.11 ±0.03 ±0.02 


1 


051 ±0 


009 


1 


00- 


-2 


00 




119 




9 





41 ±0.06 ±0.03 ±0.01 ±0.02 


1 


035 ±0 


Oil 



A recent measurement in Z — ^ t^t decays which has a larger statistical uncertainty [3] 
can also be combined with the electron and muon channels, yielding 

(t(pp -> Z -> tr ) = 75.4 ± 0.8 (stat.) ± 1.7 (syst.) ± 2.6 (lumi.) ± 0.4 (FSR) pb. 

The results may be compared with theoretical calculations similar to those used in the 
interpretation of the Z ^ analysis |2j. These calculations are performed at NNLO 

with the program FEWZ [T7] version 2.1.1 and using the NNLO PDF sets of 
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T 



Data (Stat) i LHCb 

Data (tot) 



• MSTW08 
A NNPDF21 
□ CTEQ (CTIO) ^ 



50 60 70 80 90 100 

a(Z^e*e-) [pb] 

Figure 3: Cross-section for pp — Z — )■ e"'"e~ at y/s = 7TeV measured in LHCb, shown as 
the yellow band. The inner (darker) band represents the statistical uncertainty and the 
outer the total uncertainty. The measurement corresponds to the kinematic acceptance, 
Pt > 20 GeV/c and 2 < < 4.5 for the leptons and 60 < M < 120 GeV/c^ for the dilepton. 
The points show the various theoretical predictions with their uncertainties as described 
in the text. 



MSTW08 [18], NNPDF21 [I9] or CTEQ (CTIO NNLO) [20]. In Fig. |3| we present the 
measured cross-section and in Fig. |4|^a) the measurements of the Z rapidity distribution, 
compared in each case with the three calculations. The uncertainties in the predictions 
include the effect of varying the renormalisation and factorisation scales by factors of two 
around the nominal value, which is set to the Z mass, combined in quadrature with the 
PDF uncertainties at 68% confidence level. The data agree with expectations within the 
uncertainties. 

The differential cross-section as a function of 0* is shown in Fig. |4]^b), compared with 
the predictions of QCD to NNLO. Figure |5](a) displays the ratios of these predictions to 
the measurements. The NNLO calculations tend to overestimate the data at low 0* and 
to underestimate the data at high 0*. It is expected that the 0* distribution, like that 
of Pt, is significantly affected by multiple soft gluon emissions, which are not sufficiently 
accounted for in fixed order calculations. A QCD calculation which takes this into account 
through resummation is provided by Resbos [21j|^ Another resummed calculation |22| 
has been compared with ATLAS data [23j in the central region of rapidity, but is not 
yet available for the LHCb acceptance. Alternatively, Powheg [21] provides a framework 
whereby a NLO QCD {0{as)) calculation can be interfaced to a parton shower model such 
as Pythia which can approximate higher order effects. Comparisons with these models, 
and with the LHCb version [9] of Pythia [8] are shown in Fig. ^h). The Resbos and 
Powheg distributions are normalised to their own cross-section predictions, while the 
Pythia distribution is normalised to the cross-section measured in data. It is seen that 
Resbos gives a reasonable description of the 0* distribution. Powheg shows that the 
combination of a parton shower with the 0{as) QCD prediction significantly improves 

^The P branch of Resbos is used with grids for LHC at = 7TeV based on CTEQ6.6. 
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the description of data in the low 0* region, while in the high 0* region the data are still 
underestimated. Pythia models the data reasonably well. Overall. Resbos and Pythia 
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Figure 4: Differential cross-section for pp — )■ Z — )■ e^e^ as a function of (a) Z rapidity 
and (b) (p*. The measurements based on the ^/s = 7TeV LHCb data are shown as the 
yellow bands where the inner (darker) band represents the statistical uncertainty and the 
outer the total uncertainty. NNLO QCD predictions are shown as points with error bars 
reflecting their uncertainties as described in the text. 
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Figure 5: Ratios of various QCD calculations to data for the differential cross-section 
for pp — 7- Z — )■ e"'"e~ as a function of (p*. The measurements based on the ^/s = 7TeV 
LHCb data are shown as the yellow band centred at unity where the inner (darker) band 
represents the statistical uncertainty and the outer the total uncertainty, (a) NNLO QCD 
predictions shown as points with error bars reflecting their uncertainties as described in 
the text. Small lateral displacements of the theory points are made to improve clarity, 
(b) Ratios of the predictions of Pythia, Resbos and Powheg to the data shown as 
points, with error bars that reflect the statistical uncertainties in the predictions. For 
most points, these errors are so small that they are not visible. 



12 



6 Summary 



A measurement of the pp ^> Z — )> e+e" cross-section in pp collisions at ^/s = 7 TeV using 
0.94 fb~^ of data recorded by LHCb is presented. Although the characteristics of the LHCb 
detector prevent a sharp mass peak from being seen, a clean sample of events is identified 
with less than 5% background. Within the kinematic acceptance, px > 20GeV/c and 
2 < 77 < 4.5 for the leptons and 60 < M < 120 GeV/c^ for the dielectron, the cross-section 
is measured to be 

a(pp ^ Z ^ e+e") = 76.0 ± 0.8 (stat.) ± 2.0 (syst.) ± 2.6 (lumi.) ± 0.4 (FSR) pb. 

The cross-section is also measured in bins of the rapidity of the Z and of the angular 
variable 0*. The measurements of the rapidity distribution and of the integrated cross- 
sections are consistent with previous measurements using Z decays to and t+t" 
and show good agreement with the expectations from NNLO QCD calculations. The (f)* 
distribution, related to the Z distribution, is better modelled by calculations which 
approximately include the effects of higher orders. 
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Appendix 



Table A.l: Correlation coefficients for the differential cross-section of Z — > e+e at 7 TeV 
between bins of Z rapidity, yi- Both statistical and systematic contributions are included. 



y-L bin 


2.-2.25 


2.25-2.5 


2.5-2.75 


2.75-3. 


3.-3.25 


3.25-3.5 


3.5-3.75 


3.75-4. 


4.-4.25 


2.00-2.25 


1 


















2.25-2.50 


0.47 


1 
















2.50-2.75 


0.50 


0.70 


1 














2.75-3.00 


0.51 


0.70 


0.75 


1 












3.00-3.25 


0.50 


0.69 


0.74 


0.75 


1 










3.25-3.50 


0.45 


0.62 


0.66 


0.67 


0.66 


1 








3.50-3.75 


0.35 


0.49 


0.52 


0.52 


0.51 


0.46 


1 






3.75-4.00 


0.20 


0.27 


0.29 


0.29 


0.29 


0.26 


0.20 


1 




4.00-4.25 


0.05 


0.07 


0.08 


0.08 


0.08 


0.07 


0.06 


0.03 


1 



Tabic A. 2: Correlation coefficients for the differential cross-section of Z -> e+e at 7 TeV 
between bins of (jf . Both statistical and systematic contributions are included. 



1.-2. 



4)* bin 


0.-0.05 


0.05-0.1 


0.1-0.15 


0.15-0.2 


0.2-0.3 


0.3-0.4 


0.4-0.6 


0.6-0.8 0.8-1. 


0.00-0.05 


1 
















0.05-0.10 


0.80 


1 














0.10-0.15 


0.73 


0.67 


1 












0.15-0.20 


0.63 


0.58 


0.53 


1 










0.20-0.30 


0.62 


0.58 


0.53 


0.45 


1 








0.30-0.40 


0.51 


0.48 


0.43 


0.38 


0.38 


1 






0.40-0.60 


0.46 


0.43 


0.39 


0.34 


0.34 


0.28 


1 




0.60-0.80 


0.34 


0.31 


0.29 


0.25 


0.25 


0.20 


0.18 


1 


0.80-1.00 


0.21 


0.20 


0.18 


0.16 


0.15 


0.13 


0.11 


0.08 1 


1.00-2.00 


0.23 


0.21 


0.19 


0.17 


0.17 


0.14 


0.12 


0.09 0.06 



14 



References 



R. Thorne, A. Martin, W. Stirling, and G. Watt, Parton distributions and QCD at 



LHCb, arXiv: 0808. 1847 



LHCb collaboration, R. Aaij et ai, Inclusive W and Z production in the forward 



region at4~s = l TeV, JHEP 1206 (2012) 058 



arXiv: 1204. 1620. 



LHCb collaboration, R. Aaij et ai, A study of the Z production cross-section in pp 
collisions at y/s = 7 TeV using tau final states, jarXiv: 1210 . 6289 

DO collaboration, V. M. Abazov et ai, Precise study of the Z/j* boson transverse 



momentum distribution in pp collisions using a novel technique, |Ptiys. Rev. Lett. 106 
(2011) 12200l| |arXiv: 1010.0262 



A. Banfi et ai, Optimisation of variables for studying dilepton transverse momentum 



distributions at hadron colliders, Eur. Phys. J. C71 (2011) 1600 arXiv: 1009 . 1580 



LHCb collaboration, A. A. Alves Jr. et al, The LHCb detector at the LHC, JINST 
3 (2008) S08005. 



R. Aaij et ai, The LHCb trigger and its performance, arXiv: 1211.3055 



T. Sjostrand, S. Mrenna, and P. Skands, PYTHIA 6.4 physics and manual, JHEP 



05 (200 6) 026[ |arXiv : hep- ph/0603175[ 



I. Belyaev et ai. Handling of the generation of primary events in Gauss, the LHCb 
simulation framework. Nuclear Science Symposium Conference Record (NSS/MIC)] 



IEEE (2010) 1155 



P. M. Nadolsky et ai. Implications of CTEQ global analysis for collider observables, 
[Phys. RevT D78 (2008) 013004HarXiv:0802.0007| 

GEANT4 collaboration, J. Allison et ai, Geant4 developments and applications. 



IEEE Trans. Nucl. Sci. 53 (2006) 270; GEANT4 collaboration, S. Agostinelli et al. 



GEANT4: A simulation toolkit, Nucl. lustrum. Meth. A506 (2003) 250| 



M. Clemencic et ai. The LHCb Simulation Application, Causs: Design, Evolution 
and Experience, IIJ. of Phys: Conf. Ser. 331 (2011) 032023| 

LHCb collaboration, R. Aaij et al.. Absolute luminosity measurements with the LHCb 



detector at the LHC, JINST 7 (2012) POlOlO, arXiv: 1110.2866, 



P. Golonka and Z. Was, PHOTOS Monte Carlo: a precision tool for QED corrections 
in Z and W decays, Eur. Phys. J. C45 (2006) 97, ,arXiv : hep-ph/0506026, , 



15 



[15] C. Carloni Calame, G. Montagna, O. Nicrosini, and A. Vicini, Precision electroweak 
calculation of the production of a high transverse-momentum lepton pair at hadron 



colliders, |JHEP 10 (2007) 109, arXiv: 0710. 1722 



[16] Z. Was, Private communication, http://annapurna.ifj .edu.pl/~wasm/phNLO.htm 



[17] R. Gavin, Y. Li, F. Petriello, and S. Quackenbush, FEWZ 2.0: A code for hadronic 
Z production at next-to-next-to-leading order, Comput. Phys. Commun. 182 (2011)] 



,2388, arXiv: 1011.3540, 



[18] A. Martin, W. Stirling, R. Thorne, and G. Watt, Parton distributions for the LHC, 
|Eur. Ph^^. C63 (2009) 189narXiv: 0901 . 0002, 



[19] R. D. Ball et ai, A first unbiased global NLO determination of parton distributions 
and their uncertainties, jNucl. Phys. B8 38 (2 010) 136^ |arXiv: 1002.4407| 



[20] H.-L. Lai et ai. New parton distributions for collider physics, Phys. Rev. D82 (2010) 



10740241 larXiv : 1007 . 2241t P. Nadolsky et ai. Progress in CTEQ-TEA PDF analysis, 
larXiv: 1206.33211 

[21] G. Ladinsky and C. Yuan, The Nonperturbative regime in QCD resummation 
for gauge boson production at hadron colliders, Phys. Rev. D50 (1994) 4239| 
arXiv : hep-ph/ 93 1 1341 ; C. Balazs and C. Yuan, Soft gluon effects on lepton 
pairs at hadron colliders, Phys. Rev. D56 (1997) 5558, arXiv:hep-ph/9704258; 
F. Landry, R. Brock, P. M. Nadolsky, and C. Yuan, Tevatron Run-1 Z boson data 



and Collins- Soper-Sterman resummation formalism, Phys. Rev. D67 (2003) 073016 



arXiv : h ep-ph/0212159[ 



[22] A. Banfi, M. Dasgupta, S. Marzani, and L. Tomlinson, Predictions for Drell-Yan q 
and Qt observables at the LHC, (Phys. Lett. B715 (2012) 1521 [arXrv :i205 . 4760[ 



[23] ATLAS collaboration, G. Aad et ai. Measurement of angular correlations in Drell- 
Yan lepton pairs to probe Z/'j* boson transverse momentum at y/s =7 TeV with the 



ATLAS detector, (arXiv: 1211 . 6899 



[24] S. Alioli, P. Nason, C. Oleari, and E. Re, NLO vector-boson production matched with 
shower m POWHEC, JHEP 0807 (2008y060l |arXiv : 0805 . 4802| S. Alioh, P. Nason, 



C. Oleari, and E. Re, Vector boson plus one jet production in POWHEC, [JHEP 01 
poll) 095 



arXiv: 1009.5594 



16 



